Two phase transitions in {dx2_y2 + zs)-wave superconductors 
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We study numerically the temperature dependencies of specific heat, susceptibility, penetration 
depth, and thermal conductivity of a coupled (d^2_y2 + is)-wave Bardeen-Cooper-Schriefi^er super- 
conductor in the presence of a weak s-wave component (1) on square lattice and (2) on a lattice 
with orthorhombic distorsion. As the temperature is lowered past the critical temperature Tc, a less 
ordered superconducting phase is created in d^2_y2 wave, which changes to a more ordered phase 
in {d^2_y2 + is) wave at Td. This manifests in two second-order phase transitions. The two phase 
transitions are identified by two jumps in specific heat at Tc and Td- The temperature dependen- 
cies of the superconducting observables exhibit a change from power-law to exponential behavior as 
temperature is lowered below Td and confirm the new phase transition. 
PACS number(s): 74.20.Fg, 74.62.-c, 74.25.Bt 
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The unconventional high- Tc superconductors with 
a high critical temperature Tc have a complicated lat- 
tice structure with extended and/or mixed symmetry for 
the order parameter It is generally accepted that 

for many of these high- Tc materials, the order parameter 
exhibits anisotropic behavior. However, it is difficult to 
establish the detailed nature of anisotropy, which could 
be of an extended s-wave, a pure d-wave, or a mixed 
(s -|- exp(i0)(i)-wave type. Some of the high- Tc materials 
have singlet c?-wave Cooper pairs and the order parame- 
ter has (1^2 _y2 symmetry in two dimensions H. Recent 
measurements 0] of the penetration depth A(T) and su- 
perconducting specific heat at different temperatures T 
and related theoretical analysis ||J^ also support this 
point of view. In some cases there is the signature of an 
extended s- or c?-wave symmetry |^,||. 

The possibility of a mixed (s — d)-wave symmetry 
was suggested sometime ago by Ruckenstein et al. and 
Kotliar [Q. Several different types of measurements 
which are sensitive to the phase of the order parameter 
indicate a significant mixing of s-wave component with 
a predominant (1^2 ^y2 state at lower temperatures be- 
low Tel. For temperatures between Td and the critical 
temperature Tc only the d^2_y2 state survives. There 
are experimental evidences based on Josephson super- 
current for tunneling between a conventional s-wave su- 
perconductor (Pb) and single crystals of YBa2Cu307 
(YBCO) that YBCO has mixed d ± s ov d ± is sym- 
metry B at lower temperatures. A similar conclusion 
may also be obtained based on the results of angle- 
resolved photoemission spectroscopy experiment by Ma 
et al. in which a temperature dependent gap anisotropy 
in Bi2Sr2CaCu208+a: was detected Q. The measured 
gaps along both high-symmetry directions are non-zero 
at low temperatures and their ratio was strongly tem- 
perature dependent. This observation is also difficult to 
reconcile employing a pure s- or d-wave order parameter 
and suggests that a mixed [d^2_y2 -|-exp(i0)s] symmetry 



is applicable at low temperatures. However, at higher 
temperatures one could have a pure dx2_y2 symmetry of 
the order parameter. Recently, this idea has been ex- 
plored to explain the NMR data in the superconductor 
YBCO and the Josephson critical current observed in 
YBCO-SNS and YBCO-Pb junctions (l0|. Recently, a 
new class of c-axis Josephson tunneling experiments are 
reported by Kouznetsov et al. in which a conven- 
tional superconductor (Pb) was deposited across a single 
twin boundary of a YBCO crystal. In that case measure- 
ments of critical current as a function of the magnitude 
and angle of a magnetic field applied in the plane of the 
junction provides a direct evidence of an order parame- 
ter of mixed [dx^-y^ -|-exp(i0)s] symmetry in YBCO. The 
microwave complex conductivity measurement in the su- 
perconducting state of high quality YBa2Cu3 07_5 single 
crystals measured at 10 GHz using a high-Q Nb cavity 
also strongly suggests a multicomponent superconduct- 
ing order parameter in YBCO [|2|. 

More recently, Krishana et al. reported a phase 
transition in the high- Tc superconductor Bi2Sr2CaCu208 
induced by a magnetic field from a study of the thermal 
conductivity as a function of temperature and applied 
field. Possible interpretation of this measurement could 
be the induction of a minor s or dxy component with 
a dj;2 „2 symmetry with the application of a weak field 

There have also been certain recent theoretical stud- 
ies using mixed s- and d-wave symmetries |l^,0 and it 
was noted that it is more likely to realize a stable mixed 
d -\- is state than a d -I- s state considering different cou- 
plings and lattice symmetries. As noted by Liu et al. 
]l6[ , a stable d -|- s solution can not be realized on square 
lattice. However, in the presence of orthorhombic dis- 
tortion, such a solution can be obtained. The special 
interest in these investigations was on the temperature 
dependence of the order parameter of the mixed d -I- is 
state within the Bardeen-Cooper-Schrieffer (BCS) model 
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|^,|9|. The study by Liu et al. [|6| of the order parame- 
ter on mixed d + is state also explored the effect of a van 
Hove singularity in the density of states on the solutions 
of the BCS equation. Laughlin has provided a theoreti- 
cal explanation of the observation by Krishana et al. ||l^] 
that at low temperatures and for weak magnetic field 
a time-reversal symmetry breaking state of mixed sym- 
metry is induced in Bi2Sr2CaCu208. From a study of 
vortex in a d-wave superconductor using a self-consistent 
Bogoliubov-de Gennes formalism, Franz and Tesanovic 
p5[ also predicted the possibility of the creation of a su- 
perconducting state of mixed symmetry. Although, the 
creation of the mixed state in this case is speculative, 
they conclude that a dramatic change should be observed 
in the superconducting observables when a pure super- 
conducting dx2_y2 state undergoes a phase transition to 
a mixed-symmetry state. In many cases there are evi- 
dences of a s-wave mixture with a dx2_y2 state. In view 
of this, in this study we present an investigation of this 
phase transition on specific heat, spin susceptibility, pen- 
etration depth, and thermal conductivity. 

There is no suitable microscopic theory for high-T,, su- 
perconductors. Although it is accepted that Cooper pairs 
are formed in such materials, there is controversy about 
a proper description of the normal state and the pair- 
ing mechanism for such materials In the absence of 
a suitable microscopic theory, a phenomenological tight- 
binding model in two dimensions which incorporate the 
proper lattice symmetry within the BCS formalism |^] 
has been suggested This model has been very suc- 
cessful in describing many properties of high- Tc materials 
and is often used in the study of high- Tc compounds. We 
shall use this model in the present investigation. 

In all previous theoretical studies on mixed-symmetry 
superconductors a general behavior of the temperature 
dependence of the order paramers emerged, independent 
of lattice symmetry employed . The tight-binding 

BCS model for a mixed d + is state becomes a coupled 
set of equations in the two partial waves. The ratio of the 
strengths of the s- and d-wave interactions should lie in a 
narrow region in order to have a coexisting s and d wave 
phases in the case oi d + is symmetry. As the s-wave {d- 
wave) interaction becomes stronger, the d-wave (s-wave) 
component of the order parameter quickly reduces and 
disappears and a pure s-wave (d-wave) state emerges. 

In this work we study the temperature dependencies 
of different observables, such as specific heat, suscepti- 
bility, penetration depth, and thermal conductivity, of a 
{dj.2_y2 + is)-wave BCS superconductor with a weak s- 
wave admixture both on square lattice and on a lattice 
with orthorhombic distortion and find that it exhibits 
interesting properties. A pure s-wave (d-wave) super- 
conducting observable exhibits an exponential (power- 
law) dependence on temperature. In the case of mixed 
(dx2-y2 -|-is)-wave symmetry, we find that it is possible to 
have a crossover from an exponential to power-law depen- 
dence on temperature below the superconducting critical 
temperature and a second second-order phase transition. 



For a weaker s-wave admixture, in the present study we 
establish in the two-dimensional tight-binding model (1) 
on square lattice and (2) on a lattice with orthorhom- 
bic distortion another second-order phase transition at 
T = Tci < Tc, where the superconducting phase changes 
from a pure c?-wave state for T > Td to a mixed (d + is)- 
wave state for T < Td- The specific heat exhibits two 
jumps at the transition points T = Td and T = Tc- The 
temperature dependencies of the superconducting spe- 
cific heat, susceptibility, penetration depth and thermal 
conductivity change drastically at T = Td from power- 
law behavior (typical to d state with node(s) in the order 
parameter on the Fermi surface) for T > Td to expo- 
nential behavior (typical to s state with no nodes) for 
T < Td- The order parameter for the present (d -I- is) 
wave does not have a node on the Fermi surface for 
T < Td and it behaves like a modified s-wave one. The 
observables for the normal state are closer to the super- 
conducting I — 2 state than to those for the supercon- 
ducting I — state |6). Consequently, superconductivity 
in s wave is more pronounced than in d wave. Hence as 
temperature decreases the system passes from the normal 
state to a "less" superconducting d-wave state at T — Tc 
and then to a "more" superconducting state with dom- 
inating s-wave behavior at T = Td signaling a second 
phase transition. 

The pronounced change in the nature of the supercon- 
ducting state at r = Td becomes very apparent from a 
study of the entropy. At a particular temperature the 
entropy for the normal state is larger than that for all 
superconducting states signalling an increase in order in 
the superconducting state. In the case of the present 
{dj.2_y2 + is) state we find that as the temperature is 
lowered past Td, the entropy of the superconducting 
{dx2-y2 + is) state decreases very rapidly (not shown ex- 
plicitly in this work) indicating the appearance of a more 
ordered superconducting phase and a second phase tran- 
sition. 

We base the present study on the two-dimensional 
tight binding model which we describe below. This model 
is sufficiently general for considering mixed angular mo- 
mentum states, with or without orthorhombic distortion, 
employing nearest and second-nearest-neighbour hopping 
integrals. The effective interaction is taken to possess an 
on-site repulsion (vr) and a nearest-neighbour attraction 
{va) and can be represented as 

Vkq ^Vr- Va[cOs{ka: - Qx) + 0^ COs{ky - Qy)], (0.1) 

where P — I corresponds to a square lattice, and (3^1 
represents orthorhombic distortion. On expansion, and 
keeping only the s- and d2,2_j,2-wave components of this 
interaction we have 

Vkq = -Vb - V2{coskx - 13 cos ky) {cos qx - (3 cos qy). 

(0.2) 

Here Vq = —Vr + (H-/3^)wa/2 and V2 = Wa/2 are the cou- 
plings of effective s- and d-wave interactions, respectively. 
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As we shall consider Cooper pairing and subsequent BCS 
condensation in both s and d waves the constants Vq and 
V2 will be taken to be positive corresponding to attrac- 
tive interactions. In this case the quasiparticle dispersion 
relation is given by 

ek = — 2t[cos kx + (3 cos ky — 7 cos cos ky], (0.3) 

where t and pt are the nearest-neighbour hopping inte- 
grals along the in-plane a and b axes, respectively, and 
jt/2 is the second-nearest-neighbour hopping integral. 
The energy ek is measured with respect to the surface 
of the Fermi sea. 

We consider the weak-coupling BCS model in two di- 
mensions with {dx2-y2 + is) symmetry. At a finite T, one 
has the following BCS equation 



a — — tanh — - 



2T 



(0.4) 



with 



|AqPll/2. 



We use units fcp — 1, 



where fee is the Boltzmann constant. 



The order pa- 



rameter Aq has the following anisotropic form 
Aq -I- iA2(cos qx — Pcosq y). Using the above form of Aq 



and potential (0^), Eq. ( p.4| ) becomes the following cou- 
pled set of BCS equations 



^^^^tanh^'' 



2E„ 



2T 



(0.5) 



The superconducting to normal thermal conductivity ra- 
tio Ks{T)/Kn{T) is defined by || 



Kn{T) - 1)%(1 - /q) 



Figure 1 
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Fig. 1. The s- and d-wave parameters Aq, A2 
in Kelvin at different temperatures for (dj.2_y2 + is)- 
wave models 1(a) (square lattice: full line) and 2(b) (or- 
thorhombic distortion: dashed line) described in the text 
with different mixtures of s and d waves. 

Figure 2 



A2 = ^2 E 



A2(cosfc2. — /^cosfc,,)^ , E„ ,„ 
— ^ tanh -3. 0.6) 



wher e the cou pling is introduced through Eq. In Eqs. 
(3.5) and (li) both the interactions Vb and V2 are as- 
sumed to be energy- independent constants for |eq| < Tjj 
and zero for |eql > To, where To is a mathematical cutoff 
to ensure convergence of the integrals. It should be com- 
pared with the usual Debye cutoff for the conventional 
superconductors. 

The specific heat per particle is given by |l9] 



C{T) 



2 , , / o 1 dlAqP 



dT 



(0.7) 



where /q = 1/(1 + cxp(i?q/T)). The spin-susceptibihty 
X is defined by M 



T 



(0.8) 



where ^.m is the nuclear magneton. The penetration 
depth A is defined by |l^ 



\-'^{T) = A-2(o) 



(0.9) 
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Fig. 2. Specific heat ratio C{T)/CniTc) versus T/Tc 
for models 1(a) and 1(b) on square lattice: 1(a) (full 
line) and 1(b) (dashed line). The dashed-dotted line rep- 
resents the result for the normal state for comparison. 



We solved the coupled set of equations ( |G.5| ) and ( |0.6| ) 
numerically and calculated the gaps Aq and A2 at var- 
ious temperatures for T < T^- We have performed cal- 
culations (1) on a perfect square lattice and (2) in the 
presence of an orthorhombic distortion with Debye cut 
off ksTo = 0.02586 eV {Td = 300 K) in both cases. 
The parameters for these two cases are the following: 
(1) Square lattice - {a.) t ^ 0.2586 eV, /3 = 1, 7 = 0, 
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Vq = IM, and V2 = 0.73t, = 71 K, = 25 K; (b) 
t = 0.2586 eV, /3 = 1, 7 = 0, Vq = 1.92t, and V2 = 0.73i, 
Tc = 71 K, Tci = 51 K; (2) Orthorombic distorsion - 
(a) t = 0.2586 eV, f3 = 0.95, and 7 = 0, Vb = 2.06t, and 
V2 = 0.97t, Tc = 70 K, Tci = 25 K; (b) t = 0.2586 cV, 
P = 0.95, and j ^ 0, Vq ^ 2.2t, and V2 = 0.97t, = 
70 K, Tci = 50 K. For a very weak s-wave (d-wave) cou- 
pling the only possible solution corresponds to Aq = 
(A2 = 0). We have studied the solution only when a 
coupling is allowed between Eqs. (|0.5|) and (3.6). 
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Fig. 3. 3. Specific heat ratio C(r)/C„(rc) versus 
T/Tc for models 2(a) and 2(b) for lattice with orthorhom- 
bic distortion: for notations see Fig. 2. 
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Fig. 4. 4. Specific heat jump for different Tc for 
pure s and d waves: s wave (solid line, square lattice), 
s wave (dashed line, orthorhombic distortion), d wave 
(dashed-dotted line, square lattice), d wave (dashed- 
double-dotted line, orthorhombic distortion). 

The temperature dependencies of different A's have 
been studied in details previously jl^ . Here, for com- 



pleteness, in Fig. 1 we plot the temperature dependencies 
of different A's for the following two sets of s-d mixing 
corresponding to models 1(a) (full line) and 2(b) (dashed 
line), respectively. In the coupled (dx^^-yi ~\- is) wave as 
temperature is lowered past Tc, the parameter A2 in- 
creases up to T = Tci- With further reduction of tem- 
perature, with the appearance of the parameter Aq, the 
d-wave component begins to decrease and the parameter 
A2 is suppressed in the presence of a non-zero Aq. 

We also studied the effect of orthorhombicity in this 
problem. For fixed Vq and V2, if a small orthorhombicity 
is introduced in the model, Tc (Td) decreases (increases). 
For example, if, in the model 1(a) above, we vary the pa- 
rameter f3 from 1 to 0.95, there is no coupled solution 
involving s and d waves for /? < 0.96, where only the 
pure s-wave solution survives. The coupled solution is 
observed for (3 > 0.97 in this case. For /3 = 0.97 (0.99), 
with the parameters of model 1(a), Tc = 58 K (69 K), and 
Tci = 42 K (28 K). So the effect of introducing the or- 
thorhombicity in a model is to increase Td and decrease 
Tc. As a consequence the ratio Td/Tc increases. 

In order to substantiate the claim of the second phase 
transition at T = Td, we study the temperature depen- 
dence of specific heat in some detail. The different su- 
perconducting and normal specific heats are plotted in 
Figs. 2 and 3 for square lattice [models 1(a) and 1(b)] 
and orthorhombic distortion [models 2(a) and 2(b)], re- 
spectively. The superconducting specific heat exhibits an 
unexpected peculiar behavior. In both cases the specific 
heat exhibits two jumps — one at Tc and another at Td. 
From Eq. (0.7) and Fig. 1 we see that the temperature 



derivative of |Aq|^ has discontinuities at Tc and Td due 
to the vanishing of A2 and Aq, respectively, responsi- 
ble for the two jumps in specific heat. For pure d wave 
we find that the specific heat exhibits a power-law de- 
pendence on temperature. However, the exponent of this 
dependence varies with temperature. For small T the ex- 
ponent is approximately 2.5, and for large T (T — > Tc) it 
is nearly 2. In the (c?-|-is)-wave models, for Tc > T > Td 
the specific heat exhibits d-wave power-law behavior. For 
d-wave models Cs{Tc)/Cn{Tc) is a function of Tc and /3. 
In Figs. 2 and 3 this ratio for the d-wave case, for Tc 
= 70 K, is approximately 3 (2.5) for /3 = 1 (0.95). In 
a continuum calculation this ratio was 2 in the absence 
of a van Hove singularity [^. For T < Td, we find an 
exponential behavior in both cases. 

In Fig. 4 we study the jump AC in the specific heat at 
Tc for pure s- and d-wave superconductors as a function 
of Tc, where we plot the ratio AC/C„(Tc) versus Tc. For 
a BCS superconductor in the continuum AC/C„(Tc) = 
1.43 (1.0) for s-wave (d-wave) superconductor indepen- 
dent of Tc |§,|l^. Because of the presence of the van 
Hove singularity in the present model this ratio increases 
with Tc as can be seen in Fig. 4. For a fixed Tc, the ra- 
tio AC / Cn{Tc) is greater for square lattice (/3 = 1) than 
that for a lattice with orthorhombic distortion (/3 = 0.95) 
for both s and d waves. At Tc = 100 K, in the s-wave 
(d-wave) square lattice case this ratio could be as high 
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as 3.63 (2.92). Many high-Tc materials have produced a 
large value for this ratio. 



Figure 5 
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Fig. 5. Susceptibility ratio x{T)/x{Tc) versus T/T^ 
for pure d wave (solid line, square lattice), d wave (dashed 
line, orthorhombic distortion), pure s wave (dashed- 
dotted line, square lattice), {dx2-y2+is) model 1(a) (dot- 
ted hue), {d.j;2_y2+is) model 2(b) (dashed-double-dotted 
line). In all cases Tj, « 70 K. 

Next we study the temperature dependencies of spin 
susceptibility, penetration depth, and thermal conductiv- 
ity which we exhibit in Figs. 5 — 7 where we also plot the 
results for pure s and d waves for comparison. In Figs. 5 
— 7, we show the results for pure d-wave cases on square 
lattice and with orthorhombic distortion, (d -I- is) mod- 
els 1(a) and 2(b) mentioned above, and pure s-wave case 
on square lattice. In all cases reported in these figures 
Tc ~ 70 K. For pure d-wave case we obtained power- 
law dependencies on temperature. The exponent for this 
power-law scaling was independent of critical tempera- 
ture Tc but varied from a square lattice to that with an 
orthorhombic distortion. In case of thermal conductiv- 
ity, the exponent for square lattice (orthorhombic dis- 
tortion, (3 = 0.95) is 2.2 (1.4). For spin susceptibil- 
ity, the exponent for square lattice (orthorhombic dis- 
tortion, P = 0.95) is 2.6 (2.4). For the mixed {d + is)- 
wave case, d-wave-type power-law behavior is obtained 
for Tc > T > Tci with the same exponent as in the pure 
d-wave case. For T < Td, there is no node in the present 
order parameter on the Fermi surface and one has a typi- 
cal s-wave behavior. A passage from d- to s-type state at 
Tci represents an increase in order and hence an increase 
in superconductivity @. As temperature decreases, the 
system passes from the normal state to a d-wave state 
at r = Tc and then to a s-wave-type state at T = Td 
signaling a second phase transition. 

Until now no high-Tc material has clearly shown two 
jumps in specific heat as noted in the present study. If 
this transition to a mixed state occurs at a low tem- 



perature this jump is expected to be be small and very 
high precision experimental data will be needed for its 
confirmation. The same will be needed for confirming 
this phase transition from a measurement of Knight shift 
which should show a change from a power-law to expo- 
nential dependence on temperature. However, if a transi- 
tion from a pure d-wave state to a mixed-symmetry state 
takes place in any compound, they could be identified by 
a study of the temperature dependencies of the observ- 
ables considered in this work. 
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Fig. 6. Penetration depth ratio AA(T) = [A(T) - 
A(0)]/A(0) versus T/Tc for different models. For notation 
see Fig. 5. 



Figure 7 
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Fig. 7. Thermal conductivity ratio Ks{T)/Kn{T) 
versus T/Tc for different models. For notation see Fig. 
5. 

In conclusion, we have studied the {dx2-y2 + zs)-wave 
superconductivity employing a two-dimensional tight 
binding BCS model on square lattice and also for or- 
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thorhombic distortion and confirmed a second second- 
order phase transition at T = Td in the presence of a 
weaker s wave. We liave fcept the s- and rf-wave couplings 
in such a domain that a coupled (da.2_y2 + is)-wave solu- 
tion is allowed. As temperature is lowered past the first 
critical temperature Tc, a weaker (less ordered) supercon- 
ducting phase is created in rfj.2_,,2 wave, which changes 
to a stronger (more ordered) superconducting phase in 
((ij.2_j,2 +is) wave at Td. The (d^2_y2 +zs)-wave state is 
similar to an s-wave-type state with no node in the order 
parameter. The phase transition at is also marked 
by power-law (exponential) temperature dependencies of 
C{T),x{T), AA(T) and K{T) for T > (< T^). Fur- 
thermore, the effect of orthorhombic distortion is shown 
to increase the transition temperature between these two 
superconducting phases, thus stabilizing the mixed state. 

We thank Conselho Nacional de Desenvolvimento 
Cientffico e Tecnologico and Fundagao de Amparo a 
Pesquisa do Estado de Sao Paulo for financial support. 



[1] J. G. Bednorz and K. A. Miiller, Z. Phys. B 64 (1986) 
1898. 

[2] H. Ding, Nature 382 (1996) 51; 

D. J. Scalapino, Phys. Rep. 250 (1995) 329; 

J. P. Carbotte and C. Jiang, Phys. Rev. B 48 (1993) 

4231; 

B. G. Levi, Phys. Today 49, #5, (1996) 19. 
[3] B. E. C. Koltenbah and R. Joynt, Rep. Prog. Phys. 60 

(1997) 23. 

[4] W. Hardy et al, Phys. Rev. Lett. 70 (1993) 3999; 

K. A. Moler et al, Phys. Rev. Lett. 73 (1994) 2744; 

K. Gofron et al, Phys. Rev. Lett. 73 (1994) 3302. 
[5] M. Prohammer, A. Perez-Gonzalez, and J. P. Carbotte, 

Phys. Rev. B 47 (1993) 15152; 

J. Annett, N. Goldenfeld, and S. R. Renn, Phys. Rev. B 
43 (1991) 2778; 

N. Momono and M. Ido, Physica C 264 (1996) 311; 
M. Houssa and M. Ausloos, Physica C 265 (1996) 258. 
[6] S. K. Adhikari and A. Ghosh, Phys. Rev. B 55 (1997) 
1110; J. Phys.: Cond. Mat. 10 (1998) 135; 
A. Ghosh and S. K. Adhikari, Euro. Phys. J B 2 (1998) 
31. 

[7] A. E. Ruckcnstcin, P. J. Hirschfeld, and J. Apel, Phys. 

Rev. B 36 (1987) 857; 

G. Kothar, Phys. Rev. B 37 (1988) 3664. 
[8] A. G. Sun, D. A. Gajcwski, M. B. Maple, and R. C. 

Dynes, Phys. Rev. Lett. 72 (1994) 2267; 

A. G. Sun et al., Phys. Rev. B 52 (1995) R15731. 
[9] Jian Ma et al. Science 267 (1997) 83; 

P. Chaudhari and S. Y. Lin, Phys. Rev. Lett. 72 (1994) 

1084. 

[10] J. H. Xu, J. L. Shen, J. H. Miller, Jr., and C. S. Ting, 
Phys. Rev. Lett. 73 (1994) 2492; 

Q. P. Li, B. E. C. Koltenbah, and R. Joynt, Phys. Rev. 



B 48 (1993) 437. 
[11] K. A. Kouznetsov et al., Phys. Rev. Lett. 79 (1997) 3050. 
[12] S. Sridhar et al. Physica C 282-287 (1997) 256. 
[13] K. Krishana, N. P. Ong, Q. Li, G. D. Gu, and N. 

Koshizuka, Science 277 (1997) 83. 
[14] R. B. Laughhn, Phys. Rev. Lett. 80 (1998) 5188. 
[15] M. Franz and Z. Tcsanovic, Phys. Rev. Lett. 80 (1998) 

4763. 

[16] J.-X. Zhu, W. Kim, and C. S. Ting, Phys. Rev. B 57 

(1998) 13410; 

M. Liu, D. Y. Xing, and Z. D. Wang, Phys. Rev. B 55 
(1997) 3181; 

C. O'Donovan and J. P. Carbotte, Physica C 252 (1995) 
87; Phys. Rev. B 52 (1995) 16208. 
[17] K. A. Musaelian, J. Betouras, A. V. Chubukov, and R. 
Joynt, Phys. Rev. B 53 (1996) 3598; 
Y. Ren, J. H. Xu, and C. S. Ting, Phys. Rev. B 53 (1996) 
2249; 

A. J. Mills, H. Monien, and D. Pines, Phys. Rev. B 42 

(1990) 167; 

M. T. Beal-Monod and K. Maki, Physica C 265 (1996) 

309; Phys. Rev. B 53 (1996) 5775; 

L. J. Buchholtz, M. Palumbo, D. Rainer, and J. A. Sauls, 

J. Low Temp. Phys. 101 (1995) 1079; 

M. Matsumoto and H. Shiba, J. Phys. Soc. Japan 64 

(1995) 3384; 

E. A. Shapoval, JETP Lett. 64 (1996) 625. 

[18] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. 
Rev. 108 (1957) 1175. 

[19] M. Tinkham, Introduction to Superconductivity, (Mc- 
Graw-Hill Inc., New York, 1975). 

[20] P. Monthoux, A. V. Balatsky, and D. Pines, Phys. Rev. 
B 46 (1992) 14803. 



6 



